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摘 要 : 元 胞 自动 机 模型 是 土地 利用 变化 模拟 最 有 效 和 常用 的 方法 之 一 ,其 模拟 精度 受 元 胞 大 小 .时 间 间 隔 和 邻 


域 结构 及 大 小 等 参数 的 影响 。 


以 玛 纳 斯 河流 域 绿 洲 区 耕地 为 研究 对 象 ,基于 1975 年 .1995 年 .2005 年 .2010 年 及 


2015 年 土地 利用 数据 构建 耕地 扩张 CA 模型 ,定量 分 析 CA 模型 的 参数 敏感 性 。 结 果 表 明 :(1) CA 模型 模拟 精度 具 
有 时 空 尺度 依赖 性 ,时 间 间 隔 越 短 ,精度 越 高 ,其 中 时 间 间 隔 为 5 a 的 数量 精度 ,位置 精度 及 Kappa 精 度 分 别 比 时 间 
间隔 为 20a 的 高 13.89% 、2.63% 和 8.49%; 空 间 尺 度 上 ,元 胞 大 小 越 小 ,模拟 精度 越 高 ,元 胞 大 小 为 30 m 的 Kappa 系 
数值 (0.8671) 高 于 120 m 的 模拟 结果 (0.8259)。(2) 邻 域 结构 及 大 小 模拟 结果 表明 ,摩尔 邻 域 模拟 精度 优 于 冯 诺 依 
曼 邻 域 ;同等 邻 域 下 , 邻 域 大 小 为 7x7 模 拟 精 度 比 3x3 邻 域 模拟 精度 高 ,7x7 摩 尔 邻 域 是 模拟 耕地 的 最 优选 择 。(3 ) 


耕地 模拟 结果 表明 ,2020 年 耕地 面积 较 2015 年 相 比 减 少 149.29 km ,耕地 占 比 已 由 2015 年 的 67.40% 降 至 66.18% ， 
表现 出 因 增 长 饱和 而 逐渐 减缓 的 趋势 。 此 研究 可 为 全 面 了 解 元 胞 自动 机 模型 敏感 性 提供 一 定 的 理论 基础 。 


关键 词 : CA-Markov 模 型 ; 土地 利用 变化 ; 参数 敏感 性 ; 耕地 扩张 ; 玛 纳 斯 河流 域 


土地 利用 / 覆 被 变化 (land use and cover change, 
LUCC) 是 人 类 福 社 与 自然 环境 交互 的 纽带 ,是 全 球 
变化 的 重要 组 成 部 分 '”。 土 地 利用 变化 模拟 作为 
研究 LUCC 的 主要 内 容 之 一 ,其 形式 已 由 静态 向 时 
空 动态 模拟 转化 ”, 且 逐渐 形成 了 以 经 验 -统计 模型 、 
概念 机 理 模型 及 综合 模型 为 主流 的 模拟 模型 “ ,如 
元 胞 自动 机 (cellular automata,CA ) 模 型 "多 智能 体 
(multi-agent system, MAS) 模 型 中 和 CLUE-S 模 型 ” 
等 ,上 述 模 型 在 模拟 不 同 尺 度 、 环 境 下 的 土地 利用 
变化 中 均 表现 出 一 定 的 优势 。 其 中 ,CA 模型 是 一 种 
时 间 、 空 间 、 状 态 都 离散 ,空间 相互 作用 和 时 间 因 果 
关系 都 为 局 部 的 网 格 动力 学 模型 ,其 实现 了 “ 自 下 
而 上 "模拟 复杂 系统 时 空 演 化 的 过 程 , 现 已 成 为 当 
前 国际 上 土地 利用 模拟 的 常用 模型 "。 

元 胞 大 小 .时 间 间 隔 和 邻 域 结 构 及 大 小 等 是 CA 
模型 的 重要 参数 ,对 LUCC 模拟 精度 产生 一 定 影响 ， 
因此 ,CA 模型 的 参数 敏感 性 研究 一 直 是 土地 利用 变 
化 模拟 所 重点 关注 的 问题 "。 模 拟 过 程 中 ,元 胞 大 
小 可 以 反映 数据 空间 太 度 特征 ,元 胞 大 小 的 不 同 会 
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引起 邻 域 等 参数 的 变化 ,进而 产生 不 同 的 模拟 结 
果 ; 时 间 间 隔 可 反映 元 胞 自动 机 模型 迭代 步 长 与 现 
实 志 界 时 间 的 对 应 关系 "; 邻 域 结构 及 大 小 共同 决 
定 了 中 心 元 胞 单元 的 元 胞 范围 及 数目 “2 。 然 而 ， 
现 有 研究 多 关注 单个 参数 的 敏感 性 “2 ,忽视 了 
土地 利用 变化 模拟 受 多 因素 影响 ,不 同 的 参数 会 对 
模拟 结果 造成 相应 的 不 确定 性 ,因而 ,如 何 较为 全 
面 地 定量 分 析 多 参数 的 敏感 性 进而 优化 模拟 方案 
尚 有 待 开展 深入 研究 。 

玛 纳 斯 河流 域 属 天 山北 坡 经 济 带 核心 区 域 , 为 
典型 的 山地 -绿洲 -这 漠 生态 系统 (mountain-oasis- 
desert system, MODS)。 其 中 ,山地 区 是 绿洲 和 切 漠 
的 水 源 和 物 源 ,绿洲 区 广泛 分 布 着 大 面积 农田 , 充 
漠 区 是 绿洲 形成 与 发 育 的 背景 环境 和 基质 “”。 近 
年 来 , 随 着 流域 经 济 的 迅速 发 展 及 农业 的 深度 开 
发 ,流域 内 绿洲 区 耕地 迅速 扩张 ,天 然 绿 洲 大 幅 减 
少 人 工 绿洲 显著 增加 ,荒漠 化 现象 加 剧 及 林 草 地 
退化 等 现象 频 发 ””。 鉴 于 此 ,本 文 以 玛 纳 斯 河流 
域 绿 洲 区 耕地 为 对 象 ,定量 分 析 元 胞 大 小 时间 间 
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阳 、 邻 域 结构 及 大 小 等 参数 对 CA 模型 模拟 结果 精 
度 的 敏感 性 ,确定 区 域 CA 模 型 模拟 参数 的 最 优 组 
F ,模拟 2020 年 玛 纳 斯 河流 域 耕地 分 布 ,为 提高 基 
于 元 胞 自动 机 模型 的 区 域 土地 利用 变化 模拟 精度 
提供 一 定 的 方法 依据 。 


1 研究 区 与 研究 方法 


1.1 研究 区 概况 

玛 纳 斯 河流 域 位 于 天 山北 坡 .准噶尔 盆地 南 
i% MESH AR pp 84°42'~86°33'E, 43°50’ ~45°58'N( 
1) ,地 处 干旱 区 内 陆 , 属 典型 的 温带 大 陆 性 气候 ,年 
均 降 水 量 100~200 mm, 且 时 空 分 布 不 均匀 。 流 域 
由 南 至 北 依 次 为 山地 绿洲、 荡 漠 ,垂直 景观 带 谱 完 
整 ,东西 方向 依次 分 布 着 塔 西河 、 玛 纳 斯 河 ` 金 沟 河 
等 水 系 。 WEEE Al Fr . 玛 纳 斯 县 、 沙 湾 
县 及 克拉 玛 依 市 ` 奎屯 市 .和 布 克 赛 尔 蒙 古 自治 县 
和 乌 苏 市 部 分 区 域 。 截 止 2016 年 ,流域 内 人 口 已 达 
1.55x10° A, H EAE AME 1.37x10" 50”), 

流域 绿洲 区 是 新 老 绿 洲 聚 集 地 ,其 中 , 沙 湾 EG 
纳 斯 老 绿 洲 开 发 历史 较 长 ,以 石河 子 为 代表 的 新 的 
农垦 绿洲 插花 分 布 在 两 县 之 中 ,经 过 近 40 a 的 水 土 
开发 , 旦 伴随 先进 滴灌 技术 的 广泛 应 用 ,2015 年 玛 纳 
斯 河流 域 绿洲 耕地 面积 是 1975 年 的 1 倍 以 上 ”22 ， 


85°E 


绿洲 发 展 至 今 , 现 已 成 为 新 疆 第 四 大 灌区 ,在 新 疆 
农业 生产 中 有 举足轻重 的 地 位 *。 
1.2 研究 方法 
12.1 数据 来 源 与 处 理 数据 包括 土地 利用 数据 
集 .距离 数据 集 .其 他 要 素数 据 集 ( 含 =10 CC 积温 BE 
水 量 .DEM 及 坡度 等 )。 

(1) 土地 利用 数据 集 

采用 的 遥感 影像 数据 包括 1975 年 Landsat-MSS 
影像 ,1995 年 .2005 年 .2010 年 Landsat-TM 影像 及 
2015 年 Landsat-OLI 影 像 ,成 像 时 间 均 在 7 一 8 月 , 云 
量 少 ,质量 较 高 。 采 用 1:100000 地 形 图 对 2015 年 影 
像 进行 几何 校正 后 作为 参考 图 像 ,采用 二 次 多 项 式 
转换 方程 对 其 他 影像 进行 纠正 ,平均 位 置 误差 控制 
在 0.5 个 像 元 以 内 ,采用 双 线 性 内 插 法 将 图 像 重 采 
样 至 30 mx30 m。 结 合 该 地 区 地 形 图 与 土地 利用 特 
点 ,通过 人 机 交互 式 解 译 获得 玛 纳 斯 河流 域 绿 洲 区 
五 期 耕地 数据 。 其 中 ,1975 年 和 1995 年 解 译 结果 分 
别 采用 20 世 纪 70 年 代 1:47000 黑 白 航 片 .90 年 代 彩 
色 航 片 和 地 形 图 进行 检验 ,2005 年 .2010 年 及 2015 
年 解 译 结 果 经 野外 调查 和 Google Earth 历史 高 空间 
分 辩 率 影像 检验 ,采用 分 层 随机 抽样 方法 在 研究 区 
内 随机 选取 300 个 样本 点 进行 验证 各 期 解 译 精 度 均 
达到 85% 以 上 ,满足 土地 利用 变化 研究 需求 。 


86°E 


耕地 

= 1975 年 
== 1995 年 

村 m= 200548: 

@ 县 市 mm 2010 年 

— KŘ mm 2015 年 


图 例 


45°N 


图 1 研究 区 示意 及 近 40 a 耕地 扩张 


Fig. 1 Location of the study area and nearly 40 years of cultivated land expansion 
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(2) 距离 数据 集 

根据 研究 区 内 道路 .水 系 ` 村 镇 及 已 有 耕地 等 
数据 ,基于 ArcGIS 平 台 分 别 生成 距离 道路 水系、 村 
镇 和 已 有 耕地 的 栅 格 图 层 。 

(3) 其 他 数据 集 

S10 % 积 温 、 降 水 量 数 据 及 DEM 数 据 分 别 来 源 
于 资源 环境 数据 云 平 台 (http:/www.resdc.cn/) 与 地 
理 空间 数据 云 (http:/www.gscloud.cn/) , FLAP , 210 C 
积温 .降水 量 数据 空间 分 辩 率 为 1 km, DEM 数据 空 
间 分 辩 率 为 30 m。 上 述 数 据 中 ,距离 数据 集 和 其 他 
数据 集 用 于 耕地 模拟 时 约束 性 因子 的 构建 。 

1.2.2 CA-Markov 模 型 ”马尔 可 夫 (Markov) 模 型 具 
有 无 后 效 性 和 稳定 性 , 即 在 t+1 时 刻 的 土地 利用 类 
型 状态 仅 取决 于 :时 刻 的 状态 局 ,可 由 下 式 表 示 : 

Suen =S, XP (1) 
式 中 : Seep o S, 分 别 表 示 1+1.t 时 刻 研究 区 土地 利用 
类 型 状态 向 量 ;P 是 土地 利用 状态 转移 概率 矩阵 ,其 
表达 式 及 含义 详 见 文献 [25]。 

元 胞 自动 机 模型 是 一 种 时 间 空间 、 状 态 都 离 
散 , 空 间 相 互 作用 和 时 间 因 果 关 系 都 为 局 部 的 网 格 
动力 学 模型 ” , 现 已 被 广泛 应 用 于 土地 利用 变化 模 
拟 研究 ,可 用 下 式 表示 : 

Sern Ef (Sp N) (2) 
式 中 :5 是 元 胞 的 状态 集合 ; S,、5,,, DAE ttl 
时 刻 的 系统 状态 结果 ;NN 是 元 胞 的 邻 域 范围 ;f 是 邻 
域 范 围 内 元 胞 单元 相互 作用 的 状态 转移 规则 函数 。 

Markov 模 型 着 重 探 索 时 间 尺 度 上 LUCC 变化 特 
征 与 规律 ,难以 预测 土地 利用 的 空间 格局 变化 ;CA 
模型 关注 地 块 人群. 环境 等 微观 因子 间 的 相互 作 
用 机 制 , 但 其 主要 着 眼 于 元 胞 的 局 部 相互 作用 , 存 
在 一 定 的 局 限 性 外。CA-Markov 模 型 将 马尔 可 夫 链 
长 时 间 序 列 预测 的 能 力 与 元 胞 自动 机 模拟 复杂 系 
统 空间 变化 的 优势 相 结 合 ,可 实现 动态 .多 情景 模 
拟 LUCC 的 能 力 , 且 具有 较 好 的 模拟 精度 。 

基于 IDRISI Selva 17.0% 6 ,模拟 过 程 中 ,分 别 
提取 1975—1995 4F , 1995—2005 年 及 2005 一 2010 
年 土地 利用 转移 概率 ,选取 310 CARYL DEM 坡度、 
降水 量 、 距 离 道 路 水系 .已 有 耕地 及 村 庄 的 距离 等 
空间 距离 变量 作为 模拟 耕地 的 约束 性 因子 ,分 别 选 
取 1995 年 .2005 年 及 2010 年 水 域 和 建设 用 地 作为 
模拟 的 限制 性 因子 ,设置 不 同 的 元 胞 大 小 、 时 间 间 
隔 和 邻 域 结 构 及 大 小 分 别 模拟 2015 年 耕地 现状 ,并 


Ej 2015 年 实际 耕地 (R) 进 行 对 比分 析 ,进而 选取 CA 
模型 模拟 参数 的 最 优 组 合 对 2020 年 耕地 现状 进行 
模拟 。 
1.2.3 精度 验证 

(1) 遥感 影像 分 类 精度 可 采用 Kappa 系 数 进行 
验证 ,计算 公式 如 下 : 


Kappa = 


P -P 
P-P (3) 
式 中 :P, 为 正确 模拟 比例 ;P. 是 随机 情况 下 的 正确 模 
拟 比 例 ; 忆 为 理想 分 类 情况 下 的 正确 模拟 比例 。 
Kappa 系 数 取 值 范 围 为 [0,1], 知 Kappa=1, 则 两 幅 网 
一 致 ; 当 Kappaz0.7$ 时 ,表明 一 致 性 较 高 ; 当 Kappa< 
0.4 时 ,两 幅 图 相差 较 大 。 

(2) 模拟 结果 精度 运用 IDRISI Selva 17.0 平 台 
下 CROSSTAB 模块 及 VALIDATE 模块 对 耕地 数量 
变化 和 空间 位 置 变 化 的 模拟 结果 进行 验证 。VALI- 
DATE 模 块 为 IDRISI 软 件 提供 的 一 种 改进 的 Kappa 系 
数 验 证 方法 ,能 够 同时 兼顾 数量 和 空间 位 置 的 验证 。 

(3) 为 验证 模型 模拟 结果 的 合理 性 ,研究 基于 
2005 年 和 2010 年 土地 利用 / 履 被 数据 ,利用 Markov 
模型 得 到 2005 一 2010 年 土地 利用 / 履 被 转移 面积 矩 
阵 、 转 移 概 率 抢 阵 ,并 选取 距离 道路 的 距离 .距离 水 
系 ( 含 河流 水库 湖泊) 的 距离 .距离 已 有 耕地 的 距 
离 .距离 村 镇 的 距离 210 % 积 温 数据 .降水 量 数据 、 
DEM 及 坡度 等 8 个 影响 因素 制作 转移 适宜 性 图像 
集 , 进 而 对 2015 年 耕地 进行 模拟 。 模 型 中 ,设置 元 
胞 大 小 为 30m, 时 间 间 隔 为 Sa, 迭 代 次 数 为 5, 邻 域 
大 小 及 结构 为 IDRISI 17.0 软件 默认 的 5x5 冯 诺 依 
曼 。 据 上 述 ,得 到 2015 年 耕地 模拟 结果 ,并 与 2015 
年 实际 耕地 进行 对 比分 析 。 结 果 表 明 ,模拟 结果 ( 即 
元 胞 大 小 为 30 m 的 模拟 结果 、 时 间 间 隔 为 5 a 的 模拟 
结果 ) 具 有 较 高 的 可 信和 度 , 可 用 于 模拟 2020 年 耕地 。 


2 结果 与 分 析 


2.1 元 胞 大 小 敏感 性 分 析 

本 文 对 原始 数据 进行 重 采 样 ,得 到 元 胞 大 小 分 
别 为 30 mx30 m 60 mx60 m,90 mx90 m Fil 120 mx 
120 m 的 栅 格 数据 集 。 以 2005 年 作为 起 始 年 ,将 
2005 一 2010 年 土地 利用 变化 数据 提取 的 土地 利用 
转移 概率 作为 转换 规则 。 其 中 ,CA-Markov 模型 中 
设置 滤波 器 大 小 为 5x5, 和 迭代 次 数 均 为 5, 以 此 在 4 
种 空间 尺度 下 对 2015 年 玛 纳 斯 河流 域 绿洲 区 耕地 
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图 2 不 同 元 胞 大 小 下 模拟 结果 


Fig. 2. Simulation results under different cell sizes 


进行 模拟 ,得 到 不 同 元 胞 大 小 下 模拟 结果 (图 2)。 
由 表 1 可 以 看 出 , 随 着 元 胞 大 小 的 增加 ,数量 精 
度 基本 保持 稳定 不 变 , 位 置 精 度 则 随 着 元 胞 大 小 的 
增加 也 随 之 增加 ,30 m 的 Kappa 值 高 于 其 他 元 胞 大 
小 的 Kappa 值 。 其 中 不 同 元 胞 大 小 下 模拟 结果 数量 
精度 相差 仅 为 0.02% 左 右 , 基 本 处 于 稳定 状态 。 就 
位 置 精度 而 言 , 当 元 胞 尺度 为 120 m 时 ,其 共计 包含 
16 个 30 mx30 m 的 元 胞 ,范围 覆盖 更 大 ,因此 ,其 位 
置 精 度 也 随 之 增 大 。 元 胞 大 小 为 30 m(A) 时 的 
Kappa 系数 最 高 ,达到 0.8671, 当 元 胞 大 小 从 30 mt% 
加 至 120 m(As) 时 ,Kappa 系数 总 计 下 降 4.12%。 其 
中 ,元 胞 大 小 由 30 m 增 加 至 60 m(A2) AY, Kappa 系数 
下 降 4.22% ,呈现 明显 的 下 降 趋 势 ,但 当 元 胞 大 小 由 
60 m 增 至 90 m(A;) Æ 120 mit , Kappa 系数 略微 表 
现 上 升 趋势 ,这 是 由 于 在 进行 CA-Markov 模 拟 时 ,会 


表 1 不 同 元 胞 大 小 下 模拟 结果 精度 
Tab. 1 Accuracy of simulation results under 


different cell sizes 


精度 元 胞 大 小 /m 
30 60 90 120 
数量 精度 0.8708 0.8710 0.8710 0.8709 
位 置 精度 0.9009 0.9024 0.9034 0.9037 
Kappa 系 数 0.8671 0.8249 0.8257 0.8259 


出 现 元 胞 尺 才 极限 值 的 问题 , 即 当 元 胞 大 小 超出 某 
一 极限 值 时 ,Kappa 系 数 虽 相 对 较 大 ,但 其 实际 模拟 
结果 的 精度 却 不 理想 所 。 
2.2 时 间 间 隔 敏 感性 分 析 

将 1975—2015 年 划分 为 3 种 不 同时 间 间 隔 ,长 
期 (1975 一 1995 一 2015 年 B:=20) 、 中 期 (1995 一 
2005—2015 年 , B:=10) 和 短期 (2005 一 2010 一 2015 
年 ,Bi=5) ,提取 3 种 时 间 间 隔 下 土地 利用 转移 概率 ， 
分 别 以 1995 年 .2005 年 .2010 年 为 初始 年 ,模拟 20 
al10a 和 5a 时 间 间 隔 下 2015 年 绿洲 区 耕地 现状 
(B;、B;、B,) ,并 与 2015 年 实际 结果 (R) 进 行 比 较 。 
其 中 ,CA-Markov 模 型 中 设置 元 胞 大 小 为 30 mx30m， 
滤波 器 大 小 为 5Sx5 TERR BT BA 20.10 FS, A 
同时 间 间 隔 模拟 结果 如 图 3 所 示 。 

基于 不 同时 间 间 隔 的 模拟 结果 与 2015 年 实际 
土地 利用 现状 之 间 存 在 一 定 的 差异 ( 表 2、 图 3)。 
20 a 耕地 模拟 结果 (B;) 与 2015 年 耕地 真实 结果 (R) 
相差 最 大 ,达到 2987.06 km2;5 a 耕 地 模拟 结果 (B,) 
最 接近 于 真实 结果 ,两 者 仅 相 差 1302.93 km’; 10 a ff 
地 模拟 结果 (B,) 与 真实 值 相 差 1468.41 kms HR 
表明 , 随 着 时 间 间 隔 的 缩短 ,模拟 精度 不 断 增 加 , 模 
拟 结果 更 加 趋向 于 真实 化 。 

表 2 表 明 , 随 着 时 间 间 隔 的 增加 ,数量 精度 、 位 
置 精度 及 Kappa 系 数 均 随 之 下 降 , 时 间 间 隔 为 5 a 的 
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2015 年 实际 耕地 


me 耕地 C 非 耕 地 0 30km 


图 3 不 同时 间 间 隔 下 模拟 结果 


Fig.3 Simulation results at different time intervals 


表 2 不 同时 间 间 隔 下 模拟 结果 精度 
Tab. 2 Accuracy of simulation results at 


different time intervals 


时 间 间 隔 /a 
精度 
5 10 20 
数量 精度 0.8708 0.8529 0.7319 
位 置 精 度 0.9009 0.8897 0.8746 
Kappa 系 数 0.8671 0.8422 0.7822 


模拟 结果 优 于 10 a.20 a 模拟 结果 。 
2.3 邻 域 结构 及 大 小 敏感 性 分 析 

摩尔 邻 域 和 冯 诺 依 曼 邻 域 是 元 胞 自动 机 模拟 
土地 利用 变化 常见 的 两 种 邻 域 结构 。 其 中 ,在 标准 
摩尔 邻 域 中 ,一 个 中 心 元 胞 共有 8 个 相 邻 元 胞 ,而 在 
标准 汉 诺 依 曼 邻 域 中 ,中心 元 胞 只 有 4 个 相 邻 元 
胞 。CA-Markov 模 型 中 设置 元 胞 大 小 为 30 mx30 m, 
模拟 初始 年 为 2010 年 ,和 迭代 次 数 均 为 5 ,滤波 器 分 别 
自 定 义 选取 3x3、.5x5.7x7 摩 尔 邻 域 及 汉 诺 依 曼 邻 
域 ,依照 上 述 参数 对 2015 年 绿洲 区 耕地 进行 模拟 ， 
不 同 邻 域 及 大 小 模拟 结果 见 图 4、 图 5。 

由 不 同 大 小 下 摩尔 邻 域 和 冯 诺 依 曼 邻 域 的 模 
拟 结果 ( 表 3) 可 以 看 出 ,模拟 结果 的 数量 精度 保持 
不 变 ,这 是 由 于 分 辨 率 (30 m) 均 保持 一 致 , 且 模 拟 


(b) C 


2015 年 实际 耕地 


mm 耕地 CO 非 耕 地 0 30km 


图 4 不 同 邻 域 大 小 下 摩尔 邻 域 模 拟 结果 
Fig. 4 Moore neighborhood simulation results for different 


neighborhood sizes 


2015 年 实际 耕地 
0 30km 
| | 


mm 耕地 CO 非 耕地 
图 5 不 同 邻 域 大 小 下 汉 诺 依 曼 邻 域 模拟 结果 


Fig.3 Von Neumann neighborhood simulation results under 


different neighborhood sizes 


结果 中 耕地 的 像 元 个 数 均 为 7.23x10' 个 左右 , 像 元 
个 数 相 差 较 小 ,因而 其 数量 精度 保持 稳定 。 两 种 邻 
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表 3 不 同 邻 域 及 大 小 下 模拟 结果 精度 
Tab.3 Simulation results with different neighborhood 


and size accuracy 


i CEER 摩尔 邻 域 
精度 
3x3 5x5 7x7 3x3 5x5 7x7 


数量 精度 0.8708 0.8708 0.8708 
位 置 精度 0.8990 0.9009 0.9023 
Kappa 系 数 0.8660 0.8671 0.8679 


0.8708 0.8708 0.8708 
0.8994 0.9030 0.9042 
0.8662 0.8683 0.8690 


域 下 的 位 置 精 度 和 Kappa 系数 均 随 着 邻 域 大 小 的 增 
加 而 增加 。 相 同 邻 域 大 小 下 ,摩尔 邻 域 的 位 置 精度 
分 别 比 采用 冯 诺 依 曼 邻 域 模拟 精度 分 别 高 0.04% 、 
0.21% 和 0.19% ,摩尔 邻 域 Kappa 系 数 分 别 比 汉 诺 依 
曼 Kappa 系 数 大 0.02% .0.12% 和 0.11% , 7X7 的 摩尔 
邻 域 是 模拟 区 域 耕地 扩张 的 最 优 组 合 。 
2.4 最 优 参 数 下 的 绿洲 区 耕地 预测 

据 上 述 分 析 可 知 ,元 胞 大 小 为 30 mx30 m 时间 
间隔 为 5 a、 邻 域 结构 及 大 小 为 摩尔 邻 域 7x7 的 参数 
是 模拟 玛 纳 斯 河流 域 绿洲 区 耕地 的 最 优 参数 组 合 ， 
且 模 拟 达到 了 较 好 的 效果 ,模拟 方法 及 过 程 具有 较 
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图 例 

mm 耕地 
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1 
45°N 


86°E 


图 6 2020 年 玛 纳 斯 河流 域 绿洲 区 耕地 模拟 结果 
Fig. 6 Cultivated land in Manas River Basin oasis region 


simulation result in 2020 


高 的 可 信和 度 。 结 合 以 上 参数 ,对 2020 年 玛 纳 斯 河流 
域 绿 洲 区 耕地 进行 模拟 (图 6)。 

由 图 7 可 知 ,1975 一 2015 年 耕地 表现 为 持续 增 
加 的 趋势 , 占 比 由 研究 初期 (1975 年 ) 的 35.93% 增 至 
研究 末期 (2015 年 ) 的 67.40% ,累计 涨幅 达到 
30.47% ,年 均 增 幅 为 0.79%。 这 是 由 于 自 1975 年 以 
来 ,流域 水 土 开发 逐渐 加 剧 , 大 量 未 利用 地 被 开明 ， 
且 伴 随 着 先进 滴灌 技术 的 广泛 应 用 ,耕地 扩张 趋势 


耕地 占 绿洲 区 面积 比 /2% 
上 
o 


1975 1995 2005 2010 2015 2020 


图 7 1975 一 2020 年 玛 纳 斯 河流 域 绿洲 区 耕地 
占 绿洲 区 面积 比 


Fig.7 The area ratio of cultivated land to oasis area in the 


Ou 


oasis area of Manas River Basin from 1975 to 2020 


更 为 显著 。 由 于 加 入 了 约束 性 因子 及 限制 性 因子 ， 
模拟 目标 年 (2020 年 ) 的 耕地 占 比 略 降 至 66.18% , 呈 
现 出 因 增 长 饱和 而 逐渐 减缓 的 趋势 。 


3 讨论 


元 胞 自动 机 模型 对 空间 尺度 .时间 尺度 、 邻 域 
结构 及 大 小 等 均 具 有 敏感 性 ,多 因素 共同 作用 对 模 
型 的 精准 预测 造成 了 一 定 程度 的 影响 “22 。 

(1) 通常 情况 下 ,多 数 模拟 结果 精度 均 随 元 胞 
大 小 的 增加 而 下 降 "22 , 即 元 胞 大 小 越 精细 , 模 
拟 结 果 越 高 ;反之 ,模拟 结果 精度 则 较 低 中 。 但 也 
存在 不 同 的 研究 结果 ,如 匡 浩 等 采用 正 交 试验 的 
设计 ,系统 地 分 析 了 元 胞 自动 机 模型 各 尺度 之 间 的 
关系 以 及 尺度 组 合 对 模拟 结果 的 影响 ,结果 表明 ， 
元 胞 大 小 为 60 m 的 模拟 结果 优 于 元 胞 大 小 30 m。 
这 是 源 于 其 针对 不 同 的 元 胞 大 小 采用 了 不 同 的 转 
换 规则 、 邻 域 结构 及 大 小 ,由 于 CA 模型 模拟 精度 受 
多 因素 的 共同 制约 ,因此 ,这 可 能 是 与 本 文 结 果 有 
所 差异 的 原因 所 在 。 

(2) 时 间 间 隔 作 为 CA-Markov 模 型 模拟 土地 利 
用 变化 的 重要 参数 , 与 元 胞 大 小 、 邻 域 结构 及 大 小 
一 起 作用 于 CA-Markov 模 型 ,影响 着 模型 的 模拟 精 
度 。 本 研究 分 别 选取 20 a、10a 和 5 a 三 种 时 间 间 隔 
模拟 玛 纳 斯 河流 域 绿洲 区 耕地 现状 ,结果 表明 ,时 
间 间 隔 为 5 a 的 模拟 结果 优 于 10 a 和 20 a, 且 模拟 
精度 随 着 时 间 间 隔 的 增加 而 降低 ,这 与 Sinha 及 Sa- 
mat DE Fe Ai RB, (AE ESE ER 
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却 与 本 文 结果 不 同 , 其 选择 深圳 市 龙华 镇 作为 研究 
区 ,以 1990 年 土地 利用 驱动 力 变 量 为 初始 条 件 , 分 
别 以 5a 和 10 a 为 时 间 间 隔 模拟 1995 年 与 2000 年 ， 
结果 表明 模拟 时 间 跨 度 越 长 ,精度 越 高 。 导 致 结 
出 现 差 异 的 原因 可 能 是 政策 因素 所 致 ,龙华 镇 的 城 
市 化 过 程 主要 发 生 在 1990 一 2000 年 ,于 2000 年 基 
本 达到 稳定 ,在 1990 一 1995 年 ( 即 5a 时 间 间 隔 ) 模 
拟 过 程 中 ,模拟 的 建设 用 地 与 实际 建设 用 地 虽 在 数 
量 上 大 致 相等 ,但 在 空间 分 布 上 略 有 差异 ,导致 模 
拟 误差 较 大 ;而 在 1990 一 2000 年 ( 即 10 a 时 间 间 隔 ) 
的 模拟 过 程 中 ,由 于 城市 化 过 程 已 基本 完成 ,因此 
在 空间 分 布 上 误差 较 小 ,进而 模拟 精度 得 到 大 幅 提 
高 。 而 近 40 a 玛 纳 斯 河流 域 耕地 在 不 同时 段 内 变化 
程度 均 有 所 差异 ,其 中 ,1975 一 1995 年 期 间 ,流域 水 
土 开 发 程度 急剧 下 降 ,耕地 扩张 较为 缓慢 汪 ,导致 
20 a 模拟 结果 精度 最 低 ;1995 一 2005 年 ,流域 开始 有 
计划 地 开发 水 资源 ,农田 灌溉 以 集约 、 内 涵 挖 潜 为 
EPC JO PHS SUE ,耕地 表现 
为 持续 缓慢 扩 增 呈 ,因此 ,10 a 模拟 结果 精度 优 于 
20 a 模 拟 结果 ;2005 一 2010 年 ,流域 开荒 力度 加 大 ， 
且 随 着 “三 生 用 地 ”保护 性 政策 的 实施 ,生产 生活 用 
地 (耕地 和 建设 用 地 ) 转 入 量 快速 增加 呈 ,耕地 增 速 
加 快 ,从 而 导致 5a 模 拟 结果 优 于 其 他 。 耕 地 的 变化 
受到 经 济 .社会 人口 ` 地理 和 政策 诸多 因素 的 影 
响 ,存在 着 不 确定 性 和 复杂 性 ,因此 模拟 时 间 间 隔 
的 确定 应 考虑 区 域 发展 的 历程 和 政策 等 相关 因素 。 
(3) 邻 域 结构 及 大 小 是 CA-Markov 模型 模拟 土 
地 利用 变化 的 重要 参数 之 一 ,本 文 分 别 选 取 常 见 的 
摩尔 邻 域 和 冯 诺 依 曼 邻 域 及 不 同 的 邻 域 大 小 (3x3、 
5X5 、7x7) 定 量 分 析 CA 模 型 对 其 的 敏感 性 。 结 果 表 
明 :7x7 摩 尔 邻 域 是 模拟 土地 利用 变化 时 最 优 的 邻 
域 大 小 及 结构 。 张 亦 汉 和 Pan 等 "的 研究 结果 也 
证 实 了 摩尔 邻 域 模拟 精度 优 于 冯 诺 依 曼 邻 域 ,这 是 
由 于 在 同一 空间 尺度 下 ,摩尔 邻 域 中 角 点 (左上 AE 
下 、 右 上 、 右 下 ) 对 中 心 元 胞 的 影响 高 于 冯 诺 依 曼 邻 
域 (上 .下 左右 ), 而 Sx5 汉 诺 依 曼 邻 域 包含 摩尔 邻 
域 所 有 元 胞 ( 含 左 上 左下 AEAF) WE 
会 随 着 邻 域 角 点 的 增加 而 增加 ,相同 邻 域 结构 
下 , 随 着 邻 域 半径 的 增加 ,中 心 元 胞 的 影响 力 增 大 ， 
对 模拟 结果 精度 的 提高 起 到 了 一 定 的 促进 作用 二 。 
本 研究 尚 存在 一 定 的 不 足 :@ 邻 域 结构 仅 选取 
了 规则 邻 域 (摩尔 邻 域 及 冯 诺 依 曼 邻 域 ) ,未 考虑 不 


规则 邻 域 对 模型 模拟 精度 有 何 影响 ,国外 已 有 学 者 
针对 不 规则 邻 域 进行 探索 ” ,取得 了 较 好 的 模拟 结 
Ro © CA-Markov 模 型 的 模拟 精度 一 方面 受到 模型 
本 身 参 数 的 有 影响, 男 一 方面 来 源 于 数据 源 误差 |。 
本 文 重 在 探究 CA-Markov 模型 的 参数 敏感 性 特征 ， 
考虑 到 数据 来 源 的 一 致 性 和 完整 性 , 仅 选 取 了 
Landsat 系列 影像 作为 基础 数据 ,较为 全 面 地 对 模型 
参数 敏感 性 进行 了 定量 分 析 。 然 而 ,不 同 数据 源 的 
时 间 .空间 分 辩 率 等 差异 导致 的 模拟 误差 将 会 通过 
CA 模拟 过 程 传递 。 在 进行 土地 利用 变化 模拟 时 ,可 
通过 对 比 选取 不 同 的 遥感 数据 源 的 “最 佳 尺度 ”> ， 
进而 提高 模拟 精度 。 因 此 ,未 来 的 研究 有 必要 深入 
探讨 多 源 数据 与 模型 参数 对 土地 利用 变化 模拟 精 
度 的 影响 机 制 。@) 水 资源 是 制约 干旱 区 耕地 扩张 
的 重要 因素 ,其 包含 地 表 水 (河流 湖泊 水库 等 ) 和 
地 下 水 。 研 究 初 期 (1975 一 1990 年 ) ,流域 内 修建 了 
PEGE] KRA . 夹 河 子 等 大 、 中 .小 水 库 20 多 座 , 农 
业 灌 溉 以 地 表 水 为 主 ,灌溉 方式 主要 为 畦 灌 沟 灌 ; 
1990—1998 年 灌溉 方式 发 展 为 渠 库 井 结合 灌溉 ,水 
库 灌渠 的 共同 作用 保障 了 耕地 的 持续 灌溉 ;1998 
年 后 ,流域 高 效 节 水 灌溉 ( 膜 下 滴灌 ) 技 术 大 面积 
泛 推 广 并 应 用 ,使 得 流域 实现 了 省 水 .省 劳力 ` 省 肥 
料 ,作物 增产 达 10%~20% 左 右 吕 。 随 着 人 口 数量 的 
不 断 增 加 , 流域 经 济 得 到 较 快 发 展 , 地 表 水 已 难以 
满足 工农 业 生产 及 城市 生活 ,继而 流域 内 大 量 广泛 
布 井 并 开采 地 下 水 用 以 弥补 ,由 于 井 眼 数量 难以 统 
计 , 地 下 水 开采 量 难以 获取 ,取水 量 不 尽 相 同 ,进而 
较 难 量化 。 故 此 ,文章 将 地 表 水 作为 影响 因素 对 耕 
地 扩张 进行 模拟 分 析 。 此 外 ,土地 利用 / 歼 被 变化 存 
在 一 定 的 不 确定 性 及 复杂 性 , 需 综合 考虑 自然 、 社 
会 .经 济 及 人 文 因素 对 LUCC 产 生 的 影响 ,故而 在 今 
后 的 研究 中 ,可 针对 上 述 问 题 进行 深入 研究 ,以 期 
达到 更 高 的 模拟 精度 ,更 好 地 揭示 CA-Markov 模 型 
的 敏感 性 规律 。 


4 结论 


(1) CA 模型 模拟 精度 具有 时 空 尺 度 依赖 性 。 
就 时 间 尺 度 而 言 ,模型 模拟 的 时 间 间 隔 越 短 ,精度 
越 高 ;空间 尺度 上 ,元 胞 空间 分 辨 率 越 高 ,模型 的 模 
拟 精 度 越 高 。 基 于 5 a 时 间 间 隔 、 元 胞 大 小 为 30 m 
的 模拟 结果 更 接近 于 实际 情况 。 
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(2) 邻 域 结构 及 大 小 模拟 结果 表明 , 冯 诺 依 晶 地 利用 优化 配置 中 地 理科 学 进展 , 2019, 38(1): 26-37. [Ma 
邻 域 和 摩尔 邻 域 的 位 置 精度 和 Kappa 系 数 均 随 着 邻 Bingying, Huang Jiao, Li Shuangcheng. Optimal allocation of land 
use types in the Beijing-Tianjin- Hebei urban agglomeration based 
域 大 小 的 增加 而 上 升 。 其 他 条 件 相同 的 情况 下 , E ao 
PO RE Fass, on ecological andeconomic benefits trade-offs[J]. Progress in Geog- 
尔 邻 域 模拟 结果 的 位 置 精度 和 Kappa 系数 均 比 冯 诺 raphy, 2019, 38(1): 26-37. | 
依 曼 邻 域 略 高 ,7x7 摩 尔 邻 域 是 模拟 区 域 耕地 的 最 9] Bu, 汤 蜂 , IKRE, 等. 基于 CLUS-S 模 型 的 县 域 生产 -生活 - 
优选 择 。 生态 空间 冲突 模拟 及 特征 分 析 吊 , 生态 学 报 , 2019, 39(16): 
(3) 玛 纳 斯 河 流 域 绿洲 区 2020 年 耕 地 模 拟 结 5897-5908. [Zhao Xu, Tang Feng, Zhang Pengtao, et al. Dynamic 
果 表 明 , 相 比 2015 年 耕地 面 2020 年 耕地 面积 wanes nein ee of pa A 
ai ny HeLa lp Jay ies ecological conflicts based on CLUE-S model [J]. ActaEcologica Si- 
微 下 降 ,5 a 间 共计 减少 byan km TARR 29.86 nica, 2019, 39(16): 5897-5908. ] 
km, FS RE AE KEA 32 PWT pio) 张 亦 汉 , RER ALE, 等 城市 扩张 CA 模型 的 参数 化 感性 
势 。 耕 地 是 干旱 区 人 类 生存 的 重要 基础 ,应 努力 实 APPT]. KEELER, 2018, 22(6): 1051-1059. [Zhang Yihan, Qiao- 
现 耕 地 内 部 利用 优化 与 外 部 适度 控制 并 举 之 策 , 促 Jigang, Liu Wanhua, et al. Parameter sensitivity analysis of urban 
进 绿洲 人 地 关系 的 和 谐 发 展 。 cellular automata modellJ]. Journal of Remote Sensing, 2018, 22 
(6): 1051-1059. | 
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Expansion of cultivated land in the oasis area of the Manas River 
Basin sensitivity analysis of CA-Markov model parameters 


ZHU Lei”, XIA Xin-xin”, YANG Ai-min’*, JIN Han” 
(1. College of Grassland and Environment Sciences, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, 
China; 2. Xinjiang Key Laboratory of Soil and Plant Ecological Proceses, Urumqi 830052, Xinjiang, China) 


Abstract: The cellular automaton model is one of the most effective and common methods for land-use change 
simulation. The accuracy of the simulation is affected by parameters such as cell size, time interval, and 
neighborhood structure and size. Following the land use data of the Oasis of Manas River Basin in 1975, 1995, 
2005, 2010, and 2015, this study constructed a cultivated land expansion CA model to analyze the parameter 
sensitivity of the CA model in the method of quantitative analysis. The results showed that (1) the simulation 
accuracy of the CA model depends on a spatial and temporal scale. The shorter the time interval, the higher the 
accuracy. The quantitative accuracy, position accuracy, and Kappa accuracy of a 5- year interval are 13.89%, 
2.63%, and 8.49% higher than those of a 20-year interval. At the spatial scale, the smaller the cell size, the higher 
the simulation accuracy will be. The Kappa coefficient value(0.8671) of a cell size of 30 m is higher than the 
simulation result of 120 m (0.8259). (2) The simulation results of neighborhood structure and size showed that the 
Moore neighborhood accuracy is better than Von Neumann’ s. Under the same neighborhood, the neighborhood 
accuracy with 7X7 size is higher than that of the 3X3 size. The 7X7 Moore neighborhood is the best choice for 
simulating cultivated land. (3) The simulation results of cultivated land manifested that the cultivated land 
decreased 149.29 km’ in 2020 compared with the area in 2015. The proportion of cultivated land decreased from 
67.40% in 2015 to 66.18% in 2020, which presented a gradual slow down due to growth saturation. This study can 
help provide a theoretical basis for a comprehensive understanding of the sensitivity of the cellular automaton model. 
Keywords: CA-Markov model; land-use change; parametric sensitivity; cultivated land expansion; Manas River 


Basin 


